Driving denaturation: Nanoscale thermal transport as a probe of DNA melting 
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DNA denaturation has been a subject of intense study due to its relationship to DNA transcrip- 
tion and its fundamental importance as a nonlinear structural transition. Many aspects of this 
phenomenon, however, remain poorly understood. Existing models fit quite well with experimental 
results on the fraction of unbound base pairs versus temperature, but yield incorrect results for 
other essential quantities, such as the base pair fluctuation timescales. Here, we demonstrate that 
nanoscale thermal transport can serve as a sensitive probe of the underlying microscopic physics 
responsible for the dynamics of DNA denaturation. Specifically, we show that the heat transport 
properties of DNA are altered significantly as it denatures, and this alteration encodes detailed 
information on the dynamics of thermal fluctuations and their interaction along the strand. This 
finding allows for the discrimination between models of DNA denaturation and will help shed new 
light on the nonlinear vibrational dynamics of this important molecule. 



Besides being the "molecule of life" - or perhaps be- 
cause of it - DNA lives at a unique position where its 
double-stranded structure can unravel into two single 
strands, i.e., it denaturates or "melts", via changes in 
conditions such as temperature or ionic concentration pQ . 
Local melting of DNA due to thermal fluctuations, which 
can occur well below the denaturation temperature, is 
thought to play a major role in the formation of the tran- 
scription bubble [5] • The denaturation of DNA proceeds 
via the thermal breakage (dissociation) of base pairs and 
a nonlinear interaction between base pairs sharpens this 
transition via cooperative effects in base pair unbinding 
[3H5]. As a result, the accurate description of mechanisms 
of DNA denaturation requires one to not only correctly 
account for local thermal fluctuations but also how these 
fluctuations interact along the strand. These processes 
are typically probed indirectly through the measurement 
of the DNA melting curve - the fraction of unbound base 
pairs versus temperature at equilibrium - thus hinder- 
ing the understanding of the denaturation mechanisms. 
However, thermal fluctuations and their interaction along 
the strand are the same processes that occur in thermal 
transport, which, as we shall see, can give an indepen- 
dent and more direct assessment of the proposed DNA 
denaturation mechanisms. 

During the past decade, there has been tremendous 
progress in measuring the thermal and electronic proper- 
ties of molecular systems [SJ [7] , including a recent exper- 
iment on the heat conduction of inhomogeneous DNA- 
Gold nanocomposites |S]. However, except for few stud- 
ies relevant to double-stranded DNA [91-11 1 j . no existing 
analysis of thermal transport takes into account DNA's 
large structural fluctuations that eventually result in its 
denaturation. 



* kirill@lanl.gov 

t Current address: School of Physics and Astronomy, Tel Aviv Uni- 
versity, Israel. 
1 mpz@lanl.gov 




Figure 1: Schematic of DNA between two heat reservoirs 
that probe its structure via an energy current. 



In this paper, the non-equilibrium behavior of single- 
coordinate nonlinear models of DNA are analyzed as heat 
is driven through DNA via two thermal reservoirs, as 
shown schematically in Fig. [T] The preeminent single- 
coordinate model for thermally driven denaturation of 
DNA is the Peyrard-Bishop-Dauxois (PBD) model (3H3, 
which is believed to capture the main physical processes 
behind DNA denaturation |12j . By first deriving an ana- 
lytic expression for the thermal conductance, k, we show 
that the PBD model predicts a substantial jump in k 
upon DNA denaturation. That is, the conductance of 
DNA is predicted to actually increase, despite the fact 
that the DNA becomes disordered as it is broken into two 
single strands. This change is a delicate balance of the 
model parameters and the mechanisms that it represents. 
Indeed, we show that another related model |13j . which 
also can describe the statistical properties of DNA de- 
naturation, gives qualitatively different non-equilibrium 
behavior, predicting a drop in k. Thus, by driving DNA 
out of equilibrium, the basic physical aspects of denat- 
uration - thermal fluctuations and their coupling along 
the molecule - can be probed. 

We begin with a description of the PBD model [4]. 
This model, as well as other single-coordinate models 
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for DNA, reduces the DNA strand to a one-dimensional, 
classical lattice with each base pair represented by a sin- 
gle coordinate - the collective stretching of the hydrogen 
bonds in a Watson-Crick pair. The PBD Hamiltonian 
takes on the form 



H = J2 



+ V(y n )+W(y n ,y n - l ) 



(1) 



where y n represents the stretching coordinate within the 
n th base pair. An on-site Morse potential, V(y n ) — 
D{eT aVn — l) 2 , describes the harmonic binding of the base 
pair at small y n with the angular frequency yj2Da 2 Jm, 
and unbinding at !/„ > a with the finite dissociation en- 
ergy D. The PBD nearest-neighbor potential is given 
by 

Wiy^Vn-x) = j(l + P e- a ^ + y^)(y n ~y n ^) 2 . (2) 

The PBD model has been quite successful in character- 
izing statistical properties of DNA near the dcnaturation 
transition 14. 15. . However, since it was not derived from 
a microscopic Hamiltonian, one has to validate the model 
by comparison with experimental results. In this regard, 
there is a considerable ambiguity in the results for short 
DNA strands [TU [TMT8] as well as the timescales for the 
opening of bubbles at room temperature [161 119] . Ther- 
mal transport, as we show below, can determine some 
of the main features single-coordinate models should dis- 
play (beyond equilibrium denaturation curves) and help 
settle the above discrepancies. 

We start our analysis by examining the behavior of 
the PBD model in the low (L) and high (H) tempera- 
ture limits, where it becomes harmonic with the effective 
Hamiltonian 
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+ D ^vl + it (Vn - Vn-lY 



(3) 



with /j, = L, H. Numerical studies show that the full 
PBD model, described by the Hamiltonian in Eq. ([IJ, un- 
dergoes a sharp transition corresponding to the denatu- 
ration of DNA [5]. Above the transition temperature, T c , 
the base pairs become unbound, i.e., the on-site potential 
approaches a constant leading to Dh = 0. The nearest- 
neighbor potential becomes harmonic with Kh = K. In 
the low temperature limit, the strand will again have 
harmonic behavior with the parameters of the effective 
Hamiltonian Dl = Da 2 and Kl — K(l + p) (the lat- 
ter reflects the larger nearest-neighbor coupling of the 
low temperature strand). Even at low-temperatures, the 
nonlinearity of the potential is important, but the har- 
monic limit gives the correct physical interpretation of 
the numerical results on the full model we present later. 

We introduce the concept of the thermal conductance 
ratio, R, as the ratio of the conductance between the high 
and low temperature phases R = kh/ k l- For Ohmic 
heat reservoirs strongly coupled to the first and last site 



of the strand, it can be calculated to be (see Appendix |A| 
and Refs. O [21] for details) 



R : 



2KhDl 



(4) 



which uses the lowest order term in Kl/Dl for kl (see 
All) in Appendix [A| . Since Eq. Q is for strong 
ing, it represents an extreme value for R. For the 



PBD model with the original parameters from Ref. [4], 
the corresponding ratio is Rpbd ~ 18 [22 . That is, 
the transition from low to high T results in a substantial 
jump in k due to the weakening of the on-site confining 
potential, i.e., the base pair dissociation, as T crosses T c 
from below. The removal of the on-site potential softens 
the modes of the strand, making them more efficient in 
the transport of heat. This is more transparent if one 
examines the low temperature conductance, 



K L (X 



Kj Kj 



2 ' 



(5) 



where ojl is the on-site frequency and we have taken 
Kl/Dl to be small. Equation ^ shows that the con- 
ductance increases with decreasing frequency. 

We note that the PBD model give a novel phenomenon 
compared to other transitions. For example, in the ice- 
water transition at 273 K, the conductance decreases, 
rather than increases, with a ratio of i?H 2 ~ 0.25 [23] . 
This is mainly due to the breaking of crystal structure 
across the transition that reduces the coupling between 
molecules. Another example is single-crystal Ceo, which 
exhibits a face-centered-cubic to simple-cubic structural 
transition at 260 K [21] . The ratio for this transition can 
be estimated from the experimental data to be Rc 60 ~ 
0.75. 

The thermal conductance ratio is thus important in 
identifying the underlying mechanisms of a structural 
transition. The balance of the softening of the lattice and 
the reduction of the nearest neighbor coupling determines 
both the magnitude of the conductance change, as well 
as whether it increases or decreases. Different propos- 
als for the nonlinear forms of the model will give differ- 
ent limiting harmonic forms, changing the conductance 
behavior across the temperature range not only quan- 
titatively, but also qualitatively. For example, another 
single-coordinate model of DNA is that of Joyeux and 
Buyukdagli (JB) in Ref. [15] . It uses the same on-site po- 
tential but employs an alternative nearest-neighbor po- 
tential that, although giving a very similar melting curve, 
results in a different high temperature harmonic Hamil- 
tonian (see Appendix [A|. In the strong-coupling limit, R 
changes as 

Rpbd « 18 — > R JB « 0.017. 

That is, the PBD model gives qualitatively distinct be- 
havior from the JB model because the latter looses virtu- 
ally all of its nearest neighbor coupling at high tempera- 
ture, i.e., it looks more like the ice-to- water transition in 
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character. This qualitative conclusion does not change 
even when using different parameter sets that have been 
developed [35]. Therefore, while they both agree well 
with some denaturation experiments, the conductance 
ratio can be used to directly probe the physics contained 
within their respective interaction potentials. While tem- 
peratures much higher than the transition temperature 
are beyond the validity of the models, the opposing pre- 
dictions persist within a range of temperatures, including 
around the denaturation transition. In addition, having 
the ends of the DNA "clamped" together on each of the 
heat reservoirs - one way to measure its thermal trans- 
port properties - will help the fluctuations described by 
these single-coordinate models to persist to higher tem- 
peratures. This strongly supports using the heat conduc- 
tion as a test to determine basic aspects of DNA denatu- 
ration that need to be represented in microscopic models. 

To gain further insight into the non-equilibrium behav- 
ior of DNA, we solve for the full dynamics of the PBD 
model (Eq. with the appropriate stacking interaction) 
by numerically integrating the equations of motion of the 
DNA between two Langevin reservoirs (see Appendix [B] 
for details on numerical simulations). In Fig. [2^a), we 
plot the thermal current as a function of the average 
reservoir temperature, while keeping the temperature dif- 
ference constant at a value significantly lower than the 
width of the phase transition so that the entire strand is 
always in one phase. 

As seen Fig. [2j the ability of the DNA to conduct heat 
increases substantially after it becomes denaturated ac- 
cording to the PBD model, thus confirming our analyti- 
cal prediction. Further, the conductivity at (T) = 300 K 
is 0.18 pW/KA - smaller than that in Ref. [11]. The 
difference is due to different interactions taken into ac- 
count, which result in an absence of any signatures of 
denaturation up to 350 K in Ref. [IT] . Thus, their re- 
sults do not incorporate the effect of strong nonlineari- 
ties due to denaturation/local melting. We have also cal- 
culated the thermal current for purely harmonic strands 
with the parameters corresponding to the different phases 
(as discussed below Eq. Q). The thermal conductance 
ratio for the harmonic limits is found numerically to be 
Rpbd * 6, which is lower than our analytical result due 
to the more realistic reservoir implemented in our sim- 
ulations (see Appendix |B[ ). Figure [2jc) plots the heat 
capacity C (per base pair) as a function of temperature, 
indicating the denaturation transition and the approach 
to harmonicity (where C/fcs is always unity). 

The nonlinear model, however, gives much richer be- 
havior. While at high temperature there is convergence 
of the nonlinear model to its harmonic limit, at low tem- 
peratures a non-monoticity appears, more clearly seen in 
Fig. [2jb) where we plot the thermal current as a func- 
tion of temperature on a logarithmic scale. The non- 
monotonicity, reflected in the drop of k in the tempera- 
ture range — 100 K, is mainly due to phonon scattering 
off the anharmonic on-site potential. This signature - 
conductance versus temperature - indicates the strength 
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Figure 2: Thermal transport across the denaturation transi- 
tion, (a) Thermal current through DNA as a function of aver- 
age temperature, (T) = (Tl + Th) /2, calculated while keep- 
ing the temperature difference AT = Th — Tl — 9.3 K con- 
stant. The thermal current displays a steep rise as the tem- 
perature crosses the transition temperature T c ~ 350 K. Low 
and high temperature simulations of the harmonic strands de- 
fined by Hz, and Hh are given by blue and red lines, respec- 
tively. The numerical results on the full PBD model approach 
the harmonic limits at sufficiently low and high temperatures, 
(b) Logarithmic scale plot of the thermal current. The con- 
vergence to the harmonic limit, as well as an non-monotonic 
dependence of the thermal current on temperature at low tem- 
peratures and a sharp rise near the transition are visible. The 
latter indicates that DNA may be suitable to function as a 
molecular thermal switch, (c) Dependence of equilibrium heat 
capacity, C, (normalized per length of the strand) on temper- 
ature showing the transition at T c and the convergence to the 
harmonic value Cjks = 1. The convergence to the harmonic 
value at high temperature is much more rapid than a low 
temperature, indicating that the anharmonicity of the Morse 
potential is still important even at low temperature. 



of the nonlinearity in DNA dynamics, which we will in- 
vestigate further in a future contribution. Although the 
single-coordinate models we study here were proposed 
to study DNA behavior near the denaturation transi- 
tion, such a nonlinear signature should still appear at 
low temperature, implying that a measurement of DNA 
thermal transport could also serve to define the temper- 
ature range of validity of the PBD model. It also gives 
rise to a three orders of magnitude increase of the ther- 
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mal conductance, indicating that this characteristic of 
structural phase transitions can be useful in developing 
heattronic devices such as thermal switches and thermal 
transistors |26| . 

In the above discussion we have considered only homo- 
geneous DNA strands. It is clear that heterogeneity will 
give rise to interesting effects in the conductance. Alter- 
nating sequences of AT and GC pairs, for instance, will 
tend to suppress the low temperature conductance due 
to the mismatch between parameters of base pair bind- 
ing potentials, while the high temperature conductance 
will remain roughly the same. This will result in a sub- 
stantial increase in R [2B], which will also occur due to 
random or genomic DNA. Differences between AT and 
GC pairs gives an additional experimental test of this 
class of models |5j)] . 

To conclude, we have investigated the out-of- 
equilibrium properties of single-coordinate models of 
DNA in order to elucidate the nonlinear nature of de- 
naturation. The balance of lattice softening and reduced 
coupling in different models leads to qualitatively distinct 
predictions for the thermal conductance ratio - the ex- 
perimental measurement of which will illuminate these 
fundamental aspects of DNA and thus should shed light 
on the underlying mechanisms of DNA denaturation. 
Due to the nature of the effect, experiments will be robust 
with respect to heterogeneity, defects, and contacts, un- 
like the electrical counterpart [37]. Recent experiments 
on the heat conduction of nanotubes [58] and DNA-Gold 
nanocomposites [SJ give a potential route to developing 
the experimental setup. However, it will be necessary 
to perform the experiment in conditions of high humid- 
ity so the DNA retains its B-form [55] • Such an experi- 
ment could also be done on double- and single- stranded 
DNA separately, which would remove uncertainties due 
to the lack of solvent. Fluorescence experiments to de- 
tect local openings of the DNA strand [3DJ[3T] can further 
complement the thermal transport measurements. The 
thermal conductance ratio may also be useful in under- 
standing other structural transitions, such as nanotube 
collapse. Finally, such experiments will pave the road to 
creating molecular thermal switches and open new av- 
enues in understanding and controlling thermal flow at 
the nanoscale. 
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Appendix A: Analytic derivation of the thermal 
conductance 



The thermal conductance of a classical harmonic lat- 
tice can be found analytically. Our starting point is the 
procedure of Refs. [501151]. We consider the limiting cases 
of the single-coordinate Hamiltonian, Eq. ([!]), for a lat- 
tice of length N. One already saw that the low- (L) and 
high- (H) temperature limits can be approximated by a 
harmonic Hamiltonian of the form 



2 + ^™ + 2 



(Al) 



where /i = L, H . 

To simplify the analysis, the lattice is coupled to two 
heat reservoirs at the first and last sites, which gives the 
equations of motion 



my„ 



-2(D U + Kf^yn + K^yn-i + Vn+l) 



+ (8 n ,i + 8„,n) 



d£ A{t - t')y n {t') + rj n (t) 



(A2) 



We choose the spectrum of the dissipation to be ohmic, 
A(lo) — —ijuj, with coupling 7, and the noise to be white 
noise, {t]l/h{^)vl/_h{^')) = 4nT L /ffj5(u) + u'), with 
Tl/h the l° w arL d high reservoir temperatures. This form 
for the reservoirs will satisfy the fluctuation-dissipation 
theorem. The resulting equations of motion are 



my n 



+ (S n ,l + 5 n ,N) [-7^71 (*) + Vn{t)} ■ 

The solution for the coordinates has the form 

/>oo 

y n (t) = (1/2*) / duY-^{u)f, m {uj)e^\ 



(A3) 



(A4) 



where 7) is a vector of length N with the first and last 
components being T)l/h(u)) and the rest being zero. It 
represents the coupling of the reservoirs to the ends of 
the lattice. The N x N matrix Y — — oj 2 M~ A encodes 
the solution. Here 4> nm = 2{D^+K tl )5 n<m - K a 5 n ^ m+1 - 



Mij = mSij, and An = Ann = A{uj) and 



Anm = otherwise. 

The heat current flowing into the lattice is J = 
([J_ dt'A(t — t')yi(t')]yi(t)), where the average is over 
the noise. Setting 7 = Am, the heat current becomes 

ATX 2 m 2 



Ju 



/ dujuj 2 {(Vi N - X 2 u 2 m 2 V 2 n-i) 2 

k J-00 

+ \ 2 uj 2 m 2 (V 1 , N „ 1 +V 2 , N ) 2 }~ 1 \C hN \ 2 , (A5) 



where AT = Tu — Tr is the temperature difference of 
the reservoirs, C± n 1S the cofactor of ij. jvi ancl f n,m is 
the determinant of the submatrix of (4> — ui 2 M) from the 
n-th row (column) to the m-th row (column). It follows 
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that \C hN \ 2 = Kl N - 2 and V n . m 



elements of V° m are given by J" I 

\ U 2,N Ur > 

where 



K^ m+1 Vl m . The 

L 'i N-l \ _ q-N 



2,N-1 



T 



2(1 + D^/K^-im/K^Lu 2 -1 
1 



(A6) 



We look for the eigenvalues of T of the form exp(±zg) 
with real q because those eigenvalues correspond to prop- 
agating modes. This requirement imposes that 2cos(g) = 
2(1+ D ^ / K^) — {mjK^)uj 2 . After changing variables from 
u> to q that satisfy this constraint, the final expression (for 
an infinite lattice (N — > oo)) is 



AT 



7 



2-nm 



i-2-k 

1 ' 



dq- 



sin 



2 7 2 
in A" 



5h 



cos(g) 



(A7) 



This gives, explicitly for the low and high temperature 
thermal conductance k m = J M /AT, 



k B mK 2 
4 7 3 



+ 



27 2 ^ 
mK}, 



with 



'1 + 



4 7 2 47 2 D M 8 7 4 L^ 4 7 4 i;2 



mK, 



+ 



mK 2 

h 1 



+ 



m 2 Kf, 



m 2 Kt 



(A8) 



(A9) 



With these expressions one can explicitly find the thermal 
conductance ratio R. We have verified that for reservoirs 
contacted to a single site on each end, the thermal con- 
ductance from our numerical simulations agree with our 
analytic formula to within 10 — 15%, which we attribute 
to finite size effects in the numerical simulations. 

We can take various limiting forms of these equations. 
If we define the prefactor as and a dimensionless reser- 
voir coupling as 



7m 



7 



the expressions for the conductance become 



,2 1 o„,2^V 

K„ 



(A10) 



(AH) 



with 



(A12) 

The appropriate limiting forms for our case are the fol- 
lowing. When the high temperature harmonic limit has 
no onsite potential, then the heat conductance becomes 



1 + 2 7 | - V /T+4 7 | 



(A13) 



For a low temperature limit that has a much greater 
onsite term than the nearest neighbor coupling, i.e., 
K L /D L <c 1, the heat conductance becomes 



Dr. ' 



(A14) 



which also assumes that the dimensionless coupling to the 
reservoirs is 7 l > 1. For strong coupling to the reservoirs, 
the ratio becomes 



R, 



2KhDl 
Kl ■ 



(A15) 



This is the analytic expression we use within the article. 
The strong coupling limit gives the extreme value of R. 

Now we show that the characteristic frequency in the 
PBD model is lowered as T crosses T c from below. 
For the low temperature Hamiltonian, the correspond- 
ing equation of motion is 

my n = -{2Da 2 y n + K(l + p)[(y n -y n -i) + (y n -y n+ i)}}. 

(A16) 

From the ansatz y n = yO e iwt-ikn ^ Qne obtains the 
phonon spectrum as 



muj 2 = 2Da 2 + 2K{1 + p)[l - cos(fc)]. 



(A17) 



Thus, the frequency band of phonons is yj2Da 2 jm < 
lu < yJ[2Da 2 + AK(1 + p)]/m. For the high temperature 
Hamiltonian, the equation of motion is 

my n = -K[(y n - y n _{) + (y n - y n +i)} (A18) 

and the phonon spectrum is 

raw 2 = 2K[1 - cos(fc)]. (A19) 

The frequency band of phonons is < lo < y 'AK/m. 
The two limiting Hamiltonians are both harmonic and 
the characteristic frequency is indeed lowered, and sim- 
ilar considerations apply for other models. In the low 
temperature limit, the onsite potential stiffens the DNA 
compared to the high temperature limit, which results in 
the raising of the phonon spectrum of the low temper- 
ature limit compared to the high temperature. In the 
latter, as well, the nearest neighbor coupling drops from 
K (1 + p) to K shrinking the bandwidth. This trade-off is 
responsible for the change in thermal conductance across 
the transition. If the drop in nearest neighbor coupling 
is small, then the softening will dominate, and the heat 
conductance will increase because these softened modes 
can conduct heat more effectively. 



Appendix B: Numerical Simulation Details 

To study the dynamics of the DNA out of equilibrium 
we solve numerically the Langevin equation, which de- 
scribes the dynamics of a Hamiltonian system in the pres- 
ence of thermal baths. The Langevin equation is given 



G 



by 

dW dV 

my n = -- - r n y„ + /(t), (Bl) 

oy n oy n 

where W(y n ) and V(y n ) are the potentials described in 
Eq. The DNA strand is split into three regions, the 
two ends, each of length I, serve as the Langevin ther- 
mal reservoirs at temperatures Tl and Tjj. This means 
that the friction term T n only operates for n within the 
thermal reservoirs. The fluctuating term f(t) is Gaussian 
white noise which obeys the fluctuation-dissipation rela- 
tion (/(*)/(*')) = 2T n k B T L{H) S(t - t 1 ) for the low and 
high temperature reservoirs, respectively. 

The middle region of the length M is the free DNA 
strand, which is driven out of equilibrium by the 
Langevin reservoirs when Tj, 7^ Th- The parameters 
for the simulations are M = 60 and I — 20, and the 
PBD model parameters are D = 0.04 eV, a = 4.47 A , 
K = 0.04 eV/A 2 , to = 300 u, p = 0.5, and a = 0.358 A" 1 
[2]. The strand is homogeneous DNA except for the ran- 
dom forces and friction in the Langevin regions. The 
equations of motion are integrated with the fourth-order 
Runge-Kutta method. The temperature profile of the 
strand is evaluated by defining the local temperature 
at site n as fcfiT, 1 ™ = m(y%), where the average is 
over time. The local heat current is given by J„ = 

— (|yn 3ly ^°' !< °"' 1 ^ - The simulations are performed long 
enough to allow the system to reach its steady state. 

The simulations of the heat capacity were performed 
by connecting Langevin reservoirs of identical tempera- 



ture to every site in the chain and evaluating the average 

energy per site by E — + V(y n ) + W(y n ,y n -i)\, 

where the averaging is performed over time and all the 
sites of the chain. Then, the dependence of average en- 
ergy per site on temperature, E(T), obtained by scanning 
through a set of temperature points, was numerically dif- 
ferentiated to yield the heat capacity. 

We also performed simulations on the JB model [13] , 
which has the same onsite potential as the PBD model 
and uses the nearest neighbor potential 

W{y n ,y n _ l )=K b (y n -y n _ l f + ^{l-e- b ^-^-^). 

(B2) 

The parameters in this interaction are K b = 10- 5 eV/A , 

AH = 0.44eV, and b = 0.1 A" 2 . Thus, the low temper- 
ature limit is K[ B = 2K b + bAH, with D J L B = D L . The 

coupling K[ B = 0.044 eV/A 2 is not substantially differ- 
ent from the PBD model, which has K L = 0.06 cV/A 2 . 
However, numerical simulations indicate that in the high 
temperature limit the heat conductance converges to the 
conductance of a harmonic Hamiltonian with nearest 
neighbor interaction K b {y n — y„_i) 2 . In other words, the 
quantity (y n — y n -i) 2 is sufficiently large to suppress the 
second term in Eq. ( |B2[ ). Thus, the high temperature ef- 
fective nearest neighbor coupling is Kj^ = 2Kb, which is 
substantially smaller than the low temperature coupling, 
and also substantially smaller than the high temperature 
coupling of the PBD model. 
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